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Resolving the electromagnetic mechanism of 
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Light scattering at nanoparticles and molecules can be dramatically enhanced in the 'hot spots' 
of optical antennas, where the incident light is highly concentrated. Although this effect is widely 
applied in surface-enhanced optical sensing, spectroscopy and microscopy, the underlying 
electromagnetic mechanism of the signal enhancement is challenging to trace experimentally. 
Here we study elastically scattered light from an individual object located in the well-defined hot 
spot of single antennas, as a new approach to resolve the role of the antenna in the scattering 
process. We provide experimental evidence that the intensity elastically scattered off the 
object scales with the fourth power of the local field enhancement provided by the antenna, 
and that the underlying electromagnetic mechanism is identical to the one commonly accepted 
in surface-enhanced Raman scattering. We also measure the phase shift of the scattered light, 
which provides a novel and unambiguous fingerprint of surface-enhanced light scattering. 
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Metal nanostructures can convert incident light into nano- 
scale confined and strongly enhanced optical fields 
('hot spots'), thus acting as the optical analogue of an 
antenna. When molecules or nanoparticles (objects) are placed at 
a hot spot, the Raman scattering 1-10 , fluorescence 11 ' 12 or infrared 
absorption 13,14 , can be dramatically enhanced. A variety of surface- 
enhanced Raman scattering (SERS) experiments 15-19 indicate that 
enormous signal enhancement factors of more than 10 10 are pos- 
sible, because the Raman intensity / scales with the fourth power of 
the local field enhancement / at the metal surface 3,16,20-22 , allow- 
ing even for single-molecule spectroscopy 4,5,7,23,24 . The direct and 
quantitative verification of this power law and its underlying electro- 
magnetic scattering mechanism, however, is more than a challenge 
to tackle. The difficulties arise from the complex processes typically 
involved in surface-enhanced spectroscopies, where chemical bond- 
ing and charge transfer between object and metal nanostructures 
can occur 12,25-28 . In addition, many experiments rely on measure- 
ments of samples exhibiting heterogeneously distributed hot spots 
of random field enhancement, which forces a statistical evaluation 
of the enhancement 23 . 

Here we study the light elastically scattered from an individual 
object located in the well-defined hot spot of single antennas as a 
new approach to resolve the role of the antenna in the scattering 
process. We provide experimental evidence that the scattering off 
the object is enhanced by the presence of the antenna, and that 
the underlying electromagnetic mechanism is analogous to that in 
SERS. Furthermore, the analysis of elastic scattering offers the pos- 
sibility to access the phase shift induced by the antenna, providing a 
novel fingerprint of this scattering mechanism. 

Results 

Electromagnetic mechanism of surface-enhanced scattering. In 

Figure la we illustrate the commonly accepted electromagnetic 
scattering mechanism in SERS 3,21,28 . An object (O) is located in 
close proximity to a metal nanostructure. Acting as an antenna 
(A), the metal nanostructure enhances the incoming field £ inc of 
frequency C0i by a factor / l5 generating a local field E\ oc (hot spot) 



in its proximity: 
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The metal nanostructure also enhances the scattering off the object 
at the Raman -shifted frequency co 2 by an enhancement factor f 2 . The 
field E AOA scattered off the object via the antenna is thus given by 
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The index AOA indicates the scattering path: the incoming light 
first polarizes the antenna A, the local field of the antenna polarizes 
the object O, and finally the local field of the object acts back on 
the antenna A, which radiates as a result of this interaction. We can 
approximate fi^f 2 ~fwhen the difference between C0i and co 2 is small 
and when illumination and detection direction are the same, the 
latter being a consequence of the reciprocity theorem. In this case, 
the antenna-enhanced field scattered off the object, £aoa> scales 
with the square of the local field enhancement 29 , 
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We note that equation (3) is valid particularly in the case of surface- 
enhanced elastic scattering, where C0i is exactly equal to co 2 . A 
detailed description of the elastically scattered field £aoa is provided 
in the Methods section. 

It is important to note that for a full analysis of the scattering process 
we have to consider the complex nature of the field enhancement/: 

f=\f\e^f (4) 

Here \f\ is the magnitude of the field enhancement and A(pj the phase 
shift induced by the antenna relative to the incident field. From equa- 
tions (3) and (4) it is immediately seen that the scattered intensity 



^AOaH-EaoaI 2 (typically measured in an experiment) scales with 
the fourth power of the magnitude of the local field enhancement, 
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We also find that there is a phase shift A(p AO A between the scattered field 
E AOA = I^AOAk zA(ipAOA an d me incident field £ inc , which is given by 

A(p AOA =2A(p f (6) 

The factor of 2 in equation (6) provides a distinct signature of a sur- 
face-enhanced or an antenna- enhanced scattering process. It more 
clearly reveals the role of the antenna, which is mediating both the 
illumination of the object as well as the scattering off the object. It 
is important to note that in SERS we also have to consider chemi- 
cal enhancement effects and an additional phase shift caused by the 
Raman process. These effects are not included in equations (1-6) to 
highlight the electromagnetic mechanism of the enhancement. 

Measurement principle and experimental set-up. To access both 
the magnitude and the phase shift of/, we propose an experimen- 
tal scheme as illustrated in Figure lb. We measure the field E AOA , 
which is elastically backscattered (that is, illumination and detec- 
tion directions are the same) from an individual object located in 
the well-defined hot spot of single antennas. Analysing the elasti- 
cally scattered light ensures that the field enhancements/! and/ 2 are 
exactly equal (fi=f 2 =f), and allows for interferometric measure- 
ment of the phase. Furthermore, the object should not be chemi- 
cally bound to the antenna to exclude chemical effects. In a general 
detection procedure at frequency CO, however, it is not possible to 
distinguish between the contributions E A (co) and E A q A (co) to the 
total field E^ = E A + E AOA , where E A is the field directly scattered by 
the antenna. To extract E AOA (co) from the total field E^, we modu- 
late sinusoidally the distance d between the antenna and the object 
at a frequency Q, according to d = d 0 + Adcos(Qt)- Here Ad is the 
modulation amplitude and d 0 the average antenna- object separa- 
tion. Such a modulation will not affect the direct scattering con- 
tribution E A (co), as it does not depend on the separation distance. 
However, it will modulate the contribution E A q A (co), owing to the 
nonlinear behaviour of the field enhancement with d,f=f(d). Thus, 
by demodulating the field E^ at the frequency n£l with n>l, we 
obtain the signal E m which is proportional to the field E A q A scat- 
tered off the object (for details see the Methods section): 

E n « E A0A (7) 

To demonstrate that E n (and therefore the scattered field E A q A ) is 
proportional to/ 2 , it is necessary to measure E n as a function of/. To 
this end, we measure in the following the signal E n at the hot spot of 
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Figure 1 1 Schematics of a surface-enhanced light scattering process. 

(a) Inelastic scattering process (co^a^) from an object (0) in the presence 
of a metal nanostructure that acts as an optical antenna (A), (b) Elastic 
scattering process (a)] = (02 = g>). E\ nc ((o) denotes the incident field, E^ico) is 
the field directly radiated by the antenna, and E A q A (co) is the field radiated 
by the object via the antenna. To select /Taoa> the antenna-object distance 
d is modulated at frequency Q, and the detector signal is demodulated at 
the higher harmonic frequency nQ. 
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different well-defined antennas and relate it to the calculated local 
field enhancement/. 

We note that equation (3) (and subsequently equation (7)) is valid 
for strong field enhancements/ ^>1 and does not account for the 
vectorial nature of the fields 30 . However, a more rigorous and general 
tensorial formulation (see derivation of £ AO a in the Methods section) 
proves that equation (3) indeed applies for scattering configura- 

Detector 




Figure 2 | Set-up for verifying surface-enhanced elastic light scattering. 

The set-up comprises a scattering-type scanning near-field optical 
microscope (s-SNOM, from Neaspec GmbH), which is based on an atomic 
force microscope (AFM). The apex of a conventional Si tip acts as the 
scattering object (0). The antenna and tip are illuminated at an angle of 
about 50° relative to the surface normal with the focused beam of a CO2 
laser operating at 11.1 jam. The polarization of the illuminating beam /f inc 
is parallel to the long axis of the metal rod (s-polarized light). The field 
backscattered by the antenna (A) and the tip (0) is collected and detected 
with a pseudo-heterodyne interferometer (M denotes the reference mirror 
and BS the beamsplitter). The polarizer P ensures that s-polarized light is 
detected. To obtain a modulation of the distance between tip and antenna, 
the tip is oscillating vertically at a frequency Q-300 KHz and with an 
amplitude of 70 nm. Demodulation of the detector signal at a higher 
harmonic nQ yields the amplitude and phase signals \E n \ and Acp n . 



tions where the object of interest is located at individual hot spots of 
optical antennas, as is the case in the following experiments. 

We experimentally realize the scheme proposed in Figure lb with 
the use of a scattering-type scanning near-field optical microscope 
(s-SNOM, Fig. 2) 31 ' 32 . As schematically illustrated in Figure 3a, 
the apex of the vertically oscillating Si tip mimics the object (O), while 
the metal rod acts as an antenna (A) that converts the incident infra- 
red radiation of a wavelength X= 11.1 Jim (with polarization parallel 
to the rod) into strongly concentrated fields at the rod extremities 
(hot spots) 31 . The horizontally polarized light backscattered from 
the tip-antenna configuration is detected with a pseudo-heterodyne 
Michelson interferometer 33 (Fig. 2). Demodulation of the detector 
signal at hQ yields the amplitude \E n \ and phase shift A(p n . To iden- 
tify the hot spots, we scan the tip across the antenna and record the 
sample height, as well as amplitude \E n \ and phase A(p n as a function 
of the tip position (Fig. 3b). In the amplitude map we observe the 
strongest signals \E n \ at the rod extremities, which experimentally 
proves the existence of two hot spots on the metal rod. The phase 
map shows that the signal A(p n is nearly constant across the antenna. 

Analogous to classical radio wave antennas, the resonance and 
thus highest field enhancement for Au rods is expected at a length 
L res = A e ff/2, where A e ff is the effective wavelength 34,35 . By changing 
the antenna length L at a fixed wavelength, we can thus precisely 
control the field enhancement /at the rod extremities, which is the 
key to establish a quantitative dependence between the local field 
enhancement /and the scattered field E n . We thus study in the fol- 
lowing an antenna array where the antenna length decreases from 
the top to the bottom and from the left to the right (for antenna fab- 
rication details see Methods section). In Figure 3c we show topog- 
raphy, amplitude \E n \ and phase A(p n maps of this antenna array, 
while in Figure 4a we show zoom-in images of selected antennas. 
We observe that with increasing antenna lengths L, the amplitude 
signal increases until it reaches a maximum value at a rod length of 
about L = 3.7 Jim, consistent with the theoretical antenna resonance. 
With further increasing length, the amplitude decreases. This obser- 
vation clearly reveals the resonance behaviour of the antennas 36 . The 
simultaneously acquired phase maps reveal a continuously increas- 
ing phase with increasing rod length, providing further evidence 
that we trace a resonance 37 . 





Figure 3 | Spatial mapping of surface-enhanced light scattering, (a) A Si tip (vertically oscillating at frequency U) is scanned across a Au rod antenna (A). 
Both antenna and tip (the latter mimicking the scattering object 0) are illuminated with s-polarized light. The s-polarized light backscattered from the tip- 
antenna configuration is detected interferometrically (not shown). Demodulation of the detector signal at a higher harmonic nQ yields the amplitude \E n \ 
and the phase A% at each position of the tip. (b) Topography, amplitude |£" n | and phase shift A(p n maps of a Au rod antenna. The scale bar denotes 1 Jim. 
(c) From left to right, topography, amplitude |£" n | and phase shift A(p n maps of Au rods. The length of the rods decreases from the top to the bottom and 
from the left to the right). The illumination wavelength was 11.1 jam. Demodulation was done with n = 4. The scale bar denotes 5 urn 
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O^H I Max jt ■ ^^m-n Rod length, L (urn) Phase shift, Acp f (L) (rad) 



Figure 4 | Verification of antenna-mediated elastic light scattering, (a) Amplitude \E n \ and phase shift A% maps of Au rods of different lengths. All the 
maps were recorded at an illumination wavelength of 11.1 urn Demodulation was done with n = 4. The scale bar denotes 2 urn (b) Measured intensity 
/ n °c|£ n | 2 (red symbols) and calculated field enhancement f (black solid line) at the hot spot (marked by a cross in Fig. 3b) as a function of the antenna length 
L Both quantities are normalized to their maximum value. The red solid line shows the function f 408 (L). (c) Parametric representation of log[/ n (0] and 
log[f(0] (red symbols). A linear least-square fitting of the data points (red solid line) yields a slope of 4.08. (d) Measured phase shift A(p n (red symbols) and 
calculated phase shift A(pf (black solid line) at the hot spot (marked by a cross in Fig. 3b) as a function of antenna length L. The red solid line shows the function 
1.85 Aq)f (L). (e) Parametric representation of A(p n (L) and A(pf(L) (red symbols). Linear least-square fitting of the data yields the relation Acp n = '\.85A(pf. 



To relate the scattered field E n with the local field enhancement/, 
we measure amplitude \E n \ and phase A(p n at the hot spot (marked 
by a cross in Fig. 3b) as a function of the antenna length L and com- 
pare them to the numerically calculated magnitude \f\ and phase 
Acpf of the local field enhancement. The calculations of the magni- 
tude and phase at the hot spots were performed with the use of a 
finite -difference time-domain (FDTD) method (Lumerical FDTD 
Solutions). The same wavelength and geometrical antenna param- 
eters as in the experiment were used. The dielectric values for Au 
and the refractive index of the CaF 2 substrate, n CaF2 =1.26 + 0.005/, 
were taken from Palik 38 . 

Relating scattered intensity and local field enhancement. Figure 4b 
displays the measured intensity l n °^\E n \ 2 (red symbols) and the cal- 
culated field enhancement /(black solid line) at the hot spots indi- 
cated by the cross in Figure 3b, both normalized to their maximum 
value. We clearly see the resonance behaviour, with the resonance 
appearing at L = 3.7 Jim for both I n and/. In Figure 4c we show a 
parametric representation oflog[/ n (L)] andlog[f(L)] (red symbols). 
We find a linear behaviour, revealing that the scattered intensity 
I n and the field enhancement / are related by a power law. Linear 
least-square fitting of the data points yields a slope of 4.08. By plot- 
ting / 408 (L) in Figure 4b (red solid line), we find a good agree- 
ment with the measured intensity I n (L). Our results clearly verify 
that the intensity scattered from a nanoscale object located in a hot 
spot scales with the fourth power of the local field enhancement, 
/<x/ 4 . The slight quantitative discrepancy between the experimental 
and theoretical exponent can be attributed to uncertainties in the 
antenna fabrication process and dielectric values, which have a par- 
ticularly strong impact in establishing a large -exponent power law. 

Figure 4d displays the measured phase shift A(p n (red symbols) 
and the calculated phase shift A(pf (black solid line) as a function 
of the antenna length L. As expected, both of them increase con- 
tinuously with the antenna length L (ref. 37). In Figure 4e we show 
a parametric representation of A(p n (L) and A(pj(L) (red symbols). 
Linear least-square fitting of the data yields a slope of 1.85. By plot- 
ting A(pj(L) multiplied by the factor 1.85 in Figure 4d (red solid 

4 



line), we find an excellent agreement with the measured phase shift 
A(p n (L). This factor (yielding the difference between the measured 
phase shift and the phase shift induced by the antenna) is close to 
the expected value of 2 (as predicted in equation (6)), thus clearly 
proving the double role of the antenna, which is both illuminating 
the object as well as scattering off the object. 

Antenna-enhanced light scattering is particularly efficient when 
the objects of interest are located in nanoscale antenna gaps, where 
the local field enhancement reaches maximum values 8 . To study the 
scattering mechanism in this case, we recorded the field scattered 
from the Si tip in the proximity of infrared gap antennas (Fig. 5a). 
The amplitude and phase maps of the resonant antenna are shown 
in Figure 5b and c, clearly visualizing the hot spot in the antenna 
gap. Owing to the capacitive coupling across the gap 31 , this hot spot 
exhibits a much stronger field enhancement and field confinement 
compared with that observed in single rods (Fig. 3b). From a zoom- 
in image (Fig. 5b), we measure for the hot spot a full width at half 
maximum of about 50 nm. Following an identical procedure as the 
one applied in Figure 4, we establish the relation between the field 
scattered from the tip and the local field enhancement exactly at the 
gap centre. To this end, we mapped gap antennas with a fixed gap and 
varying arm length L, and measured the field \E n \ and the phase shift 
A(p n at the gap centre as a function of L. We also calculated the mag- 
nitude |/| and phase Acpyof the field enhancement at the gap centre. 
From a parametric representation of \og[I n (L)] andlog[f(L)] (Fig.5d), 
we again find that the intensity scattered by the tip scales with nearly 
the fourth power of the field enhancement. The difference between 
the exponent of 4.56 determined from our study and the expected 
value of 4 can be assigned to uncertainties in the antenna fabrication 
process, yielding variations in the antenna dimensions (for example, 
length, width and gap size) with respect to the design values used for 
the near-field calculations. Particularly difficult to reproduce by elec- 
tron-beam lithography on the insulating CaF 2 substrate is the gap, 
which is precisely where the local field enhancement in Figure 5d 
is evaluated. In Figure 5e we show a parametric plot of A(p n (L) 
and A(pj(L), yielding clear evidence that the phase shift of the scat- 
tered field £aoa is twice the phase shift of the local field at the hot 
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In conclusion, our results provide a rigorous verification of the 
widely accepted electromagnetic enhancement mechanism in sur- 
face- and antenna- enhanced spectroscopy, more than 30 years after 
it was proposed to explain the electromagnetic contribution in SERS. 
We also establish that the enhancement mechanism for elastic and 
inelastic (Raman) scattering is the same, and thus generalize the elec- 
tromagnetic mechanism of SERS. Furthermore, the verification of the 
phase shift involved in the elastic scattering process introduces a new 
paradigm to characterize and manipulate surface-enhanced coherent 
scattering processes, which might be of fundamental and technologi- 
cal importance in optical information and sensing applications. 
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Figure 5 | Elastic light scattering in the hot spot of infrared gap antennas. 

(a) A Si tip (vertically oscillating at frequency Q) is scanned across a 
Au rod antenna (A). Both antenna and tip (the latter mimicking the 
scattering object O) are illuminated with s-polarized light. The s-polarized 
light backscattered from the tip-antenna configuration is detected 
interferometrically (not shown). Demodulation of the detector signal 
at nQ yields the amplitude \E n \ and phase shift A% at each position of 
the tip. (b) Zoom-in amplitude \E n \ and phase shift A% images of the 
gap region in a resonant gap antenna. A full width at half maximum of the 
hot spot of about 50 nm is measured. The scale bar denotes 100 nm. 

(c) Topography, amplitude \E n \ and phase shift A% maps of a resonant gap 
antenna. The maps were recorded at an excitation wavelength of 11.1 jam 
and demodulation was done with n = 4. The scale bar denotes 2jim. 

(d) Parametric representation of log[/ n (0] and log[f(0] (red symbols). 

A linear least-square fitting of the data points (red solid line) yields a slope 
of 4.56. (e) Parametric representation of A(p n (L) and Acpf(L) (red symbols). 
Linear least-square fitting of the data yields the relation A<p n = 1.99A(pf. 

spot (A(p n (L) = 1.99A(pj(L)). These observations confirm the results 
obtained with single rods (Fig. 4) and thus corroborate the electro- 
magnetic mechanism in antenna-enhanced light scattering. 

Discussion 

Our findings might have particular implications in surface- 
enhanced infrared absorption (SEIRA) spectroscopy, where objects 
adsorbed on rough metal films or antennas exhibit an enhanced 
infrared fingerprint spectrum. SEIRA — actually yielding extinction 
spectra— is typically assigned to an enhanced infrared absorption 
of the objects. The enhanced infrared scattering from the objects is 
barely considered. However, the scattering from the object might be 
a significant contribution to the extinction, as it scales with/ 4 , while 
the infrared absorption scales only with f 2 . Future interpretation 
of SEIRA thus might need to consider also the scattering process. 
On the other hand, the strongly enhanced scattering cross-section 
might be applied for the development of surface-enhanced infrared 
scattering spectroscopy. 



Methods 

Derivation of Eaoa- I n the following we show that the elastic light scattering from 
an object at the hot spot of an optical antenna, £aoa> scales with the fourth power 
of the local field enhancement. We consider a general elastic light scattering proc- 
ess similar to the one shown in Figure lb, where incoming and scattered radiations 
occur at frequency co. Both antenna (A) and object (O) are characterized by means 
of their corresponding polarizabilities a A and (Xq, with a A ^>(XQ. When an incident 
field E inc (co) impinges on the antenna and object, a multiple scattering process 
occurs. The resulting total scattered field Ej(co) can be expressed as a coherent sum 
of the contribution of the different scattering orders as: 

(8) 



ontnoution or me airrerent scattering oraers as: 
E r (co) = E A (co) + E Q (co) + E AO (co) + E OA (co) + E 



I AO A (») + ■■ 



From now on, an implicit co dependence will be assumed in all of the expressions. 
In equation (8) the term 

E A = G A a A E inc (9) 

is the field radiated directly by the antenna with G A = G A (r^><=°) being the cor- 
responding Greens function evaluated in the far field, for example, at a distant 
detector. The term 

E 0 =G 0 a 0 E inc (10) 

is the field radiated directly by the object, where G Q = Go(r^°°) is the Greens 
function describing the scattering from the object. The term 



E OA ~ G O a oAr 



(ID 



is the field radiated directly by the object after being illuminated by the local field 
^loc = fE]nc produced by the antenna. The term 



E AO ~ G kf a O E inc 
is the field scattered from the object via the antenna. The term 
£ AOA = G A/ a oAnc 



(12) 



(13) 



describes the field scattered from the object via the antenna, after being illuminated 
by the local field produced by the antenna. The direct scattering off the object Eq 
in equation (8) can be neglected because the polarizability of the object c/q is, by 
definition, smaller than the polarizability of the antenna a A . Thus, Eq is negligible 
when compared with E A . We also note that the exact total radiation in equation 
(8) should include all the successive interaction events in the scattering process. 
Nevertheless, the small value of (Xq allows for a truncation of Ej at the first order 
of Oq, that is, after the term £aoa> as shown in equation (8). Moreover, if the field 
enhancement is large (/^>1), as it is the case in a hot spot of an antenna, the con- 
tributions Eq A and £ao (which scale linearly with/) to Ej are much smaller than 
E AOA (which scales with the square off). The total radiated field Ej (equation (8)) 
can be thus reduced to: 



'-E, +E 



AOA 



(14) 



The field £aoa i n equation (13), which is the one studied in this communication, 
can be evaluated considering the tensorial nature of the polarizability a Q (diagonal 
tensor of components (Xqh with i = x, y, z) and the field enhancement f (tensor 
of components^- with i,j = x, y, z). With an incident field E- mc x polarized in the 
x-direction (denned here to be parallel to the long axis of the antenna, as in our 
experimental situation) we obtain: 

^0,xxfxxfxx E inc,x ~*~ ^0,yyfyxfyx E inc,x ~*~ 
0 



c AOA,a 



0 



(15) 



where fy is the field enhancement factor of the /-component induced by the 
j- component of the incident field E- mc j. 
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In the case of an isotropic polarizability of the object, cXq xx = cxo^y = 0Cq zz , equa- 
tion (15) yields the well-known scalar power law E A q Ax °c [(f xx ) 2 + (f yx ) 2 + (f zx ) 2 ] =f 2 
as described by equation (13). In case of an anisotropic object (as the tip in our 
experiment), the general solution of E AOAx (equation (15)) prevails. However, 
at the extremity of single dipolar antennas and at the gap between two coupled 
antennas (the two situations considered in this communication), one component 
of the enhancement dominates over the others, fxx^fyx'fzx i n me g a P °f a coupled 
antenna andf zx ^>f xx ,f yx on top of the single antenna extremities. In these cases, 
we can consider f~f xx ovf~f zx , respectively, and E AOA x can be expressed as 

E AOA,x (®) « G Afxx a O,xxfxx E mc,x 00 f ( 16 ) 

at the antenna gap, or 

E AOA,x " G Afzx a O,zzfzx E inc,x 00 f 2 ( 17 ) 

on the extremity of an isolated antenna. Equations (16) and are identical to the 
scalar description (equation (13)) of surface -enhanced light scattering. 

Mathematical relationship between £„ and Eaoa- m tne following we analyse 
the dependence of the n-th order demodulated detector signal E n (co) on the field- 
enhancement factor. E n (co) is obtained as the n-th order coefficient of the Fourier 
transform: 

E n {co) oc \e t (co) e- inQt dt = \(E A (G)) + E AOA (cD))e- Mt dt 

= jE AOA (co)e- Mt dt (18) 

Substituting E A q A (co) according to equation (13), we obtain: 

E„(w) 0= jG A a 0 E inc (<o) ■ f(d) ■ f(d) ■ e~ innt At 
= G A a 0 E lac (w)jf(d) ■ f(d) ■ e- inQt dt 
Normalizing the integral in equation (19) by/=/ 0 =/(d 0 ), we find: 

E n (co) = G A a 0 E inc m /o 2 J(/(rf)// 0 ) 2 e^'dt W 
By denning the normalized n-th Fourier coefficient % n (do) as 

Ud 0 ) = j(f(d)/fo) 2 e- innt dt (21) 
we can express E n (co) in terms of/ 0 and E, n (do) as 

E n (co) = G A a Q E inc (co)fZZ n (do ) = E AOA (co)^ n (d 0 ) (22) 

In Supplementary Figure SI we show that for linear dipole antennas t; n (d 0 ) is 
independent of the antenna length. Thus, in equation (22), ^ n (d 0 ) can be taken as a 
constant for every single antenna at a particular distance d 0 , yielding 

E n (co)<xE AOA (co)<xf 2 (23) 

Antenna fabrication. We fabricated Au rods with a nominal length L between 2 and 
5 urn and a fixed width of 100 nm on a double- side-polished CaF 2 (100) substrate by 
high-resolution electron-beam lithography. A bilayer polymethyl methacrylate (with 
495,000 molecular weight at the bottom and 950,000 molecular weight at the top) 
was used as electron- sensitive polymer. The Au film with a thickness of 40 nm was 
thermally evaporated in high vacuum on top of a 3-nm thick Ti layer, electron-beam 
evaporated for adhesion purposes. A lift-off process yields the Au nanorods. 
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